# Compare District 2 with 3 other districts on dead tree metrics
comparison_districts <- c(2, 5, 19, 26) # District 2 + 3 comparison districts
district_comparison <- TREE_DIST |>
st_drop_geometry() |>
filter(CounDist %in% comparison_districts) |>
group_by(CounDist, borough) |>
summarise(
total_trees = n(),
dead_trees = sum(tpcondition == "Dead", na.rm = TRUE),
stumps = sum(tpcondition == "Stump", na.rm = TRUE),
dead_and_stumps = dead_trees + stumps,
pct_dead = (dead_and_stumps / total_trees) * 100,
current_maples = sum(grepl("maple", genusspecies, ignore.case = TRUE), na.rm = TRUE),
pct_maples = (current_maples / total_trees) * 100,
.groups = "drop"
) |>
arrange(desc(pct_dead))
# Display comparison table
kable(district_comparison,
digits = 2,
col.names = c("District", "Borough", "Total Trees", "Dead", "Stumps",
"Total Dead+Stumps", "% Dead", "Current Maples", "% Maples"),
caption = "Tree Health Comparison: District 2 vs. Selected Districts")